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Gause's principle of competition between two species is studied when 
one of them is sterile. We study the condition for total extinction in the 
niche, namely, when the sterile population exterminates the native one by an 
optimal use of resources. A mathematical Lotka-Volterra non linear model 
of interaction between a native and sterile species is proposed. The condition 
for total extinction is related to the initial number M of sterile individuals 
released in the niche. In fact, the existence of a critical sterile-population 
value M c is conjectured from numerical analysis and an analytical estimation 
is found. When spatial diffusion (migration) is considered a critical size 
territory is found and, for small territory, total extinction exist in any case. 
This work is motived by the extermination agriculture problem of fruit flies 
in our region. 
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In ecological systems Gause's exclusion principle is widely accepted [1-5]. 
Originally it was deducted from competition between Paramecium cauda- 
tum and Paramecium aurelia [1-3]. Nevertheless, it applies to many other 
situations. For instance, in reference [6] was conjectured that Neanderthal 
extinction in Europe was consequence of Gause's principle. In a formal point 
of view, it states that two competing species cannot coexist in the same eco- 
logical niche. In this framework, it is assumed that the strong species fulls 
completely the niche and the weak disappears (exclusion). We remark that 
this principle is limited in the sense that it applies when re-adaptation, mi- 
gration, or genetic changes does not exist. This principle seems very intuitive 
in natural environment or for species in laboratories; but what is the situ- 
ation with genetically prepared sterile populations ?. To be more explicit, 
consider the well known problem related to extermination of native fruit flies 
by genetically sterilized fruit flies [7,8]. The two species exist in the same 
ecological niche when the sterile population is released. Before the interac- 
tion, we assume that the native species fulls the niche in a stable way. In 
some geographic regions and for optimal initial conditions, native fruit flies 
can be exterminated by the sterile population. Namely, in this case both 
species disappears and the principle must be reformulated as 

Principle of Gause : Two competing species cannot coexist in the same 
ecological niche and at least one of the species disappear. 

Namely, it contains explicitly the possibility of total extinction (both 
species). This formulation of the principle includes all strategy of extermi- 
nation with genetically altered species [7,8]. 

In this paper we consider a mathematical non linear model of competition 
between a native species with number of individuals N(t), and other sterile 



2 



with number of individuals given by M(t). Explicitly, we are interested on 
the mathematical conditions for total extinction in the ecological niche. The 
structure of the article is the follows: first we present a predator-predator 
non linear model for the variable N and M (equations 1-3), and its stability 
analysis. Numerical simulation confirms the stability analysis and the exis- 
tence of a threshold M c where total extinction exist. We give an analytical 
estimation of this threshold value (equation (5)). Near to the critical value, 
the behavior of the extinction time r ext is studied. This extinction time is, 
after our numerical calculation, related to an critical exponent v (equation 
(6)). Finally, the case including diffusion is considered, here we found the 
existence of a critical size territory L c were total extinction holds (for any 
initial condition of the species M). Some possible generalizations are stressed 
at the end of the paper. 



To be explicit, consider the Lotka-Volterra type evolution equations with 
interaction 

^£ = - a 'M - 5NM, (1) 
dt w 

dN 

^— = NF(N) - 5NM, (2) 

where a' is the death-rate of the sterile population, 8 is the interaction param- 
eter, and the function F(N) describes the population growing of the native 
species when interaction does not exist. For instance, when F(N) = a — N 
(a is a constant) we obtain the usual Verlhust, or logistic, equation. Remark 
that the stability of the point (N = 0, M = 0) depends of the sign of F(0). 
In fact, when F(0) < this point becomes stable and the possibility of total 
extinction exists in accordance with the Gause's principle. Moreover, this 
condition of stability of (0, 0) seems reasonable if we think that species needs 
a minimal social structure, or genetical diversity, to survival (i.e. a minimal 
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number of individuals). We stress that the dynamical systems (1-3) is irre- 
versible, for instance, a Lyapounov function L associated to the systems is 
just L = M(t). 



To carry-out explicit our calculations we consider the model where 

F(N) = (a-N)(N-(3), (3) 

with a and (3 positives parameters (a < (3). Since -F(O) < 0, the point 
(N = 0, M = 0) is stable and total extinction would be expected. The linear 
analysis of (1-3) shows that the point (N = a, M = 0) is unstable (saddle) 
and (N = (3, M = 0) is a stable focus. Figure 1 shows the stability diagram 
for our equations. So, the sterile population M disappears and, depending 
on the initial conditions, total extinction would exist in the niche. Namely, 
the systems has two atractors, the first (0,0) related to total extinction and 
the other (/3, 0) related to survival of species N. 



Figure 1: A sketch of the critical points associated with equations (1-2) 
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The explicit question which we are concerned in this paper is the follows: 
if initially the native species number is iV = f3 (the stable point without 
interaction) then, after released M sterile individuals, when have we total 
extinction?. Namely, before the interaction, the native species is in the niche 
in a stable way. After, M a sterile individuals are released and the interaction 
process produces a competitive struggle. Here we ask about the minimal 
population M of sterile individuals producing total extinction in the niche. 
In fact, if the sterile population is not enough then they death and total 
extermination do not occur. 



Numerical calculations confirm the existence of a critical value M c and 
when M Q > M c total extinction exist in the niche. Figure 2 shows the time 
behaviors of N(t) for different initial value M a , of the sterile species released 
in the niche. There is a critical value for the initial condition M Q related to 
total extinction. A criterion for total extinction is depending on the initial 
number of sterile population M a and given by 

M > 2.7 ^"^ . (4) 

This criterion is established as follows: from (1) and (3), we have that M(t) = 
exp (-«' -6<N> t )t, where < N > t = (1/t) f N(r)dr . When t -> oo, 
and assuming total extinction, we expect an exponential decaying behavior 
for M. So, an important fraction of the decaying process, assumed slow, 
is reached when M ~ M e -1 . If now we stressed that not other stationary 
point (excepting (0, 0) exist in (2), we obtain the criterion (4). We have used 
the maximum value of the function F(N) given by (3). We remark that 
this is a coarse criterion, nevertheless, it works in accord with our numerical 
simulations. For instance, the figure 2 describes extinction when M > 270 
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in accordance with (4). This is true also for other parameter values. 



Figure 2: A numerical simulations of equations (1-3) 

The criteria (4) can be generalized easily to a system with arbitrary dis- 
tribution F(N) in (1-2). Namely, by imposing the inequality F max < <5M D e _1 
(with e ~ 2,7). 



The figure 3 shows a simulation of the final native population N(t = oo) 
for different initial condition M Q of the sterile population. Clearly there is 
a critical value M c which separates the survival and extermination regime. 
From equation (4), a first approximation for this critical values is 

c ~ 7 — AS — • ^ 
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Figure 3: The final distribution of native population M 

Moreover, figure 2 suggests that near to this critical value the extinction 
time T ext depends on (M Q — M c ). This is so because when M Q — > M c + the 
extinction time must be infinity. Explicitly we expect a behavior like 

Text ~ ~, " 777; M Q > M c , (6) 

ext {M -M c ) v ' y ' 

where v is a unknown parameter. The evaluation of this critical exponent 
requires a computational work beyond of the scope of this paper. It will be 
do elsewhere. The conjecture (6) is reinforced by numerical calculation. In 
fact, using the same parameter values that in figure 2, and the definition 

r e -i = -- hm Ln(N(t)/f3), (7) 

t t~*co 

we see that the existence of the critical exponent v is confirmed numerically 
(figure 4). 
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Figure 4: The extinction time r ext 

Now we discuss briefly the incorporation of migration to the set of evolu- 
tion equations (1-2). In fact, total extermination could be also carried out by 
diffusion process. In some cases unstable points become stable by diffusion 
in limited territories. Thus in a model where (0,0) is unstable, i.e. only one 
species survive, diffusion would change this instability and total extinction 
takes place. We add the diffusion terms Dm|j and D^^^- to (1) and 
(2) respectively. Namely, consider the pair of reaction-diffusion evolution 
equations 

dN d 2 N 

— = NF(N)-5NM + D N —. (9) 

The linear analysis of stability for the stationary point (0, 0) can be do in 
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the usual way, namely, consider the small perturbation 

M = + r], 



(10) 



N = + e, (11) 

where the variables i] and e are assumed small. Equations (8) and (9), give 
the first order linear equations 

^ = -a'M, + D M p x , (12) 
de d 2 e 

_ = f(0)e+Bk _, ( 13 ) 

where we assume F(0) > 0, corresponding to the unstable case when mi- 
gration is not present. For finite territory, solutions like e ~ e^sinkx can 
be visualized. The relationship between the stability parameter uj and the 
wavenumber k is given by 

u = F(0) - k 2 D N , (14) 



and clearly for k > y F(0)/ Dm the point (0,0) becomes stable and total 
extinction in the niche exist. Since k ~ with L the territory size, equation 



(14) defines a critical size territory L c ~ y D^/ F(0) where total extinction 
holds. Namely, for all size territory L < L c total extinction exist. 



In conclusion: the principle of Gause was generalized to consider the case 
of biological struggle when one competing species is sterile. In fact, under 
appropriate conditions, total extinction could occur in the niche. Most of the 
agriculture competitive extermination method are carried-out assuming this 
principle. For instance, this is the case of extermination program of fruit flies 
with sterile flies irradiated in laboratories [9-10]. We have presented a simple 
model which has total extinction in the niche in some cases. Conjectures 
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related to a critical sterile-population (number of individuals) producing total 
extinction were pointed-out with a coarse criterion (4). This conjectures are 
based in our numerical simulation of the model. The role of migration was 
briefly discussed and the possibility of total extinction from diffusion was 
also explored for small territories. 



To ending, we note that our model can be extended to incorporates some 
modifications. Particularly we are thinking about generalizations like: 

(i) Periodic variation of coefficients. In fact, in the extinction fruit flies 
programs, daily, season or El niho (ENSO), oscillations must be considered. 

(ii) Sexual selection. Many extermination programs are based on sexual 
selection, namely, sterile-male released in a given niche. It creates interaction 
between sterile-male and fertile-female which becomes directly related to the 
evolution of the native-male. Such models require a phase-space greater than 
two. 

(iii) Many random factors are present in real niche. For instance, humid- 
ity, temperature, wind, etc. These factors can be incorporated in our model 
by introducing adequate sthocastic process for the parameter (a, (3, or 5). 
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